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Inflation is now an accepted paradigm in standard cosmology, with its predictions consistent 
with observations of the cosmic microwave background. It lacks, however, a flrm physical theory, 
with many possible theoretical origins beyond the simplest, canonical, slow-roll inflation, including 
Dirac-Born-Infeld inflation and k-inflation. We discuss how a hierarchy of Hubble flow parameters, 
extended to include the evolution of the inflationary sound speed, can be applied to compare a 
general, single field inflationary action with cosmological observational data. We show that it is 
important to calculate the precise scalar and tensor primordial power spectra by integrating the 
full flow and perturbation equations, since values of observables can deviate appreciably from those 
obtained using typical second-order Taylor expanded approximations in flow parameters. As part 
of this, we flnd that a commonly applied approximation for the tensor to scalar ratio, r ~ IGcst, 
becomes poor (deviating by as much as 50%) as Cs deviates from 1 and hence the Taylor expansion 
including next-to-leading order contribution terms involving Cs is required. By integrating the 
full flow equations, we use a Monte-Carlo-Markov-Chain approach to impose constraints on the 
parameter space of general single held inflation, and reconstruct the properties of such an underlying 
theory in light of recent cosmic microwave background and large-scale structure observations. 



I. INTRODUCTION 

Even though inflation explains our observable uni- 
verse remarkably, we have very little understanding of 
the physical mechanism responsible for the acceleration 
during inflation. Finding a physical basis for inflation 
is likely to help in understanding particle physics at 
very high energies. Many models of inflation are moti- 
vated by supergravity, the string landscape and D-branes 
[H El [31 m [5] . Therefore understanding inflation may also 
be useful in testing string theory [6, 7, 8, 9, lOj. 

Recent advances in precision cosmology provide valu- 
able constraints on the cosmological density perturba- 
tion, which is essential to understand the inflationary sce- 
nario. Improved measurements of the temperature and 
polarization anisotropies of the cosmic microwave back- 
ground (CMB) dH m m [H [151 and data from 
large scale structure surveys [T71 [TSl [THj |M1 HI] together 
characterize the primordial spectrum of fluctuations to 
fine detail. 

Using observations to constrain the primordial power 
spectrum one can reconstruct properties of the under- 
lying theory guiding the physics of the inflationary era 

[36, 37, 38, 39, 40} [ID [42l [43l [Ml [45l [46l [IB [48l [49]. 

In order to consider what observations tell us, without 
any theoretical bias, we need to reconstruct the entire 
inflaton action, instead of just the inflaton potential or 
a speciflc kinetic term, since theories of inflation, such 
as those arising from the Dirac-Born-Infeld (DBI) action 
[T] or from k-inflation [SHI [51] [511 [52j j allow the presence 
of nonminimal kinetic terms. A hierarchy of derivatives 
of the Hubble expansion factor, "flow parameters", dur- 
ing inflation was developed as a technique to reconstruct 
canonical inflation [HI [S31IS51 [55] . This was recently ex- 
tended to DBI inflation by also considering derivatives 
of the inflaton sound speed [10] and, through an addi- 
tional derivative of the Lagrangian, to a general single 



field action in [15] . 

The paper is organized as follows. In Sec. jlT] we re- 
view the background evolution equations and the fiow 
formalism in general, single field inflation. In Sec. |III| 
we discuss how we calculate the exact primordial scalar 
and tensor perturbation spectra by integrating the flow 
equations, and review the ability of approximate Taylor 
expansions about a pivot point to describe physical ob- 
servables such as the tilt, running and tensor to scalar 



IV 



ratio in the general inflationary scenario. In Sec 
we present the main findings of the paper, cosmological 
constraints on the general, single field inflationary action 
in light of current CMB temperature and polarization 
power spectra and three point temperature correlation, 
large-scale structure power spectrum and supernovae lu- 
minosity distance constraints. We consider constraints 
within the observed rang e of ph ysical scales 10"'* Mpc^^ 



< fc < 1 Mpc"* in Sec. 



IV B 



as well as general action 



reconstruction over the extended inflationary history in 
Sec. |IV C[ In Sec. |V]we draw together our flndings and 
discuss implications for the future. 



II. THE HUBBLE FLOW FORMALISM 



Consider the general lagrangian C{X, (f) of a single 
scalar fleld inflationary model. Here X = ^d^cjid^cf) is the 
canonical kinetic term. The pressure and energy density 
are given by. 



p{x,<t>) = /:(A,0), 

p{X,4,) EE 2XCx-C{X,. 



(1) 

(2) 



where Cx = dC/dX. We assume that the null energy 
condition p + p> 0, is satisfied, such that. 



(3) 
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The adiabatic sound speed for the propagation of in- 
homogeneities, c^, is defined as, 



and 
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Px 



1 



XX 



C 



X 



(4) 



We measure the extent of inflation using the variable N^, 
which denotes the number of e-folds before the end of 
inflation. We choose Ne to increase backwards in time 
from the end of inflation, i.e.. 



dN^ = -Hdt, 

a{t) ' 



In 



(5) 
(6) 



where a(t) is the scale factor at any time t, and ten is the 
time at the end of inflation. 

We can define three physical slow-roll parameters to 
describe time derivatives of the Hubble parameter and 
sound speed, 
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Hcs^ 



(7) 



where a dot represents a derivative with respect to time, 
t. Note that these parameters are independent of a scalar 
field definition. They depend upon £, and combinations 
of X and derivatives of £ with respect to X and cj) that 
are invariant under a scalar field redefinition. The accel- 
eration equation can now be written as 



a 
a 



requiring e < 1 for inflation to occur. 

The slow-roll approximation requires that 



(8) 



(9) 



where = de/dN^, etc. 

In order to describe an action beyond the slow-roll as- 
sumption, one can define an infinite hierarchy of "flow 
parameters", as used extensively for canonical inflation 
[Ml [S5J US] and extended to DBI inflation [TD], and to 
a general action in [46j. For a general action, with a 
general scalar field definition, the evolution is described 
by three hierarchies of the fiow parameters, dealing with 
derivatives with respect to the scalar field of the Hubble 
constant (77), the speed of sound (c^), and Cx- These 
parameters are in general all dependent on the explicit 
choice of and as discussed in [JB] , actions reconstructed 
using this formalism can map onto each other through a 
scalar field redefinition. In this paper we impose a spe- 
cific scalar field choice, such that Cx = c^^ , consistent 
with canonical and DBI inflation, to alleviate this degen- 
eracy. This leaves us with only two distinct hierarchies 
of flow parameters, 



(10) 
(11) 
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(13) 



for / > 1. Here a prime denotes derivative with respect 
to (f), M^i = I/SttG, and i/I'+il = d'+^ H / d<j)^+^ etc. The 
combination of parameters, 2^A — k = 2e — r/, is invariant 
under scalar field redefinition. 
Using 



dN, 



Cx H' 



(14) 



we can write the evolutionary paths of the fiow parame- 
ters as a set of coupled first order differential equations 
with respect to Nf., 



ex = -e(2e - 2^A -I- k) = -£77, 
Kx = — K(e — ^ A + 2k) + e^a, 



and for Z > 1, 



(15) 
(16) 



(17) 



ax 



= -'a[{l - l)e -{I- + + 1)k] + '+ia. (18) 



In this paper we consider two scenarios in which infla- 
tion is driven by the inflationary flow equations, one in 
which the end of inflation arises from when e = 1, and 
one in which inflation does not end on its own (e ^ 1), 
but may be brought on, for example, by the behavior of 
a second scalar fleld. 



III. PRIMORDIAL PERTURBATIONS 

In this section we discuss the generation of primordial 
power spectra in single field infiation: we summarize the 
evolution equations for the scalar and tensor perturba- 
tions in Sec. 



Ill A the choice of initial conditions in Sec. 
|IIIB| and how the exact power spectra are calculated 
through evolving the fiow equations in Sec. [Ill C[ We 
also review the approximate expressions for the power 
spectra in terms of the flow parameters in order to com- 
pare them with the exact power spectra we use for the 
analysis in Sec. |IV| 



A. Calculating the power spectrum 

The evolution of the scalar perturbations in the metric, 
ds^ = (1 + 2<i>)df - (1 - 2<i>)a^{t)-f,jdx'dx^, (19) 



3 



are typically described in terms of the Bardeen parame- 
ter, c, 



5p + 3p ^ 



2p $ 

3(p + p)^ ■ 3(p + p)i?' 

and specifically its spectral density, 



27r2 



(20) 



(21) 



while the tensor perturbations can be characterized by 
a metric with goo = ^1; zero space-time components 



We can decompose these per- 



goi = 0, and 5g^j = , 
turbations into two independent polarization modes, de- 
noted -t- and X , since gravitational waves are both trans- 
verse and traceless. Writing the Fourier modes as hk.+ 
and hk^xi the spectral density of tensor fluctuations, Vh 
can be written as, 



(22) 



The evolution of ( and h± can be calculated concisely 
through considering two alternative Mukhanov variables. 



faMpi\ 



where 



(23) 
(24) 

(25) 



To determine the full evolution of the power spectrum, 
we need to numerically integrate the mode equations in 
Uk and Vk- Written in terms of the number of e- foldings 
N„. these are 



dN? ^ ^' dN, 
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with, 
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Uk = 0, (26) 
Vk = 0, (27) 



W = 2 
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Following [57l [58], the scalar spectral density, V( is 
given by 



= 2 



2i/-3 



r(^) 



r(3/2) 



(1 



2;y-l 



2ttV2X 



%H 



(29) 



where 



3 ri K, 

^ = 2+^+2 + 2- 



(30) 



The tensor spectral density is 



r(M) 



r(3/2) 



(1 - ef^-' 



V2H 



ttM, 



pi 



k=aH 



where 



2-e 1 



(l-e)2 4 



(31) 



(32) 



B. Scalar and tensor perturbation initial conditions 

We assume the standard choice of initial conditions for 
the mode functions Uk and Vk, the Bunch-Davies vacuum. 



Uki—Cgkr oo) 
Vk{—kT — > oo) = 



1 



^ — iCs kr 



x/2^ 



2k 



(33) 
(34) 



As in [44] we note that we cannot use these condi- 
tions directly in order to solve the mode equations nu- 
merically since we cannot impose these conditions in 
the infinite past. We need to initialize the mode func- 
tions at sufliciently early times, which we choose as the 
number of e-folds before the end of inflation at which 
{csk/aH)/{\ — e — k) = 50 for scalar perturbations, and 
{k/aH)/{\ — e) = 50 for tensor perturbations. Note that 
the results are insensitive to the precise condition chosen. 

We deflne the ratio of the Hubble radius to the proper 
wavelength of fluctuations for the scalar and tensor per- 
turbations, respectively, as 



Csk 

k 



Then, 



dye 
dT 
de 

dye 
dn 

dye 



and 



dyh 

dT 

de 
dyh 



-Csk{l - e - k), 

1 de 

ye{\-e~K)dN,' 

1 dn 
veil -e-K)dNj 



-fc(l-e), 

1 de 
yh{l-e)dNl 



(35) 
(36) 

(37) 
(38) 

(39) 

(40) 
(41) 
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The initial conditions for each mode need to be set at 
early times that correspond to large y,^ and yh, when the 
scalar/tensor mode is well within the horizon. We see 
from the above equations that at large and yh, e(2/c)' 
and e(y/i) are approximately constant. Then we 
can integrate the equations in and yh to get, 



yc 

Vh 



-kT{l - e). 



(42) 
(43) 



Using this in ( 33 1 and ( 34 ) we get the initial conditions 



duk 



1 



\/2csfc 



Vi=vci 



iyci/0--<^i-Ki) 



and, 



Vk{yhi) 

dvk 
dN~. 



JVhi/il-ei) 



2k 



2k 



(44) 
(45) 

(46) 
(47) 



C. Calculating the primordial power spectrum 

We set the initial conditions for each scalar fc-mode at 
the number of e-folds, A^e, for which y^,;/(l — e^ — k^) = 50. 
We then integrate the mode equation (26 1 to find Uk{Ne) 



as we go forward in time. At each instant we obtain the 
value of z from (25 1. The flow parameters are simultane- 



ously integrated using their first order differential equa- 
tions (15)-(18l, as are the speed of sound and aH, using 



dcs 
dK 



~KCs 



d{aH) 



(1 - e)aH. (48) 



We find the scalar power spectrum (for each value of fc, 
using C from equation (23 1 in (21 1 ) as we evolve forward 
in time until the power spectrum freezes out, the condi- 
tion for which we set as [dln'P(^/dlna] < 10^^. At this 
level the accuracy of the power spectrum calculation is 
at least as good as the accuracy of the other numerical 
calculations required when obtaining the CMB and mat- 
ter power spectra predictions in CAMB, described in Sec. 
IV This allows efficient computational calculation of the 
primordial spectrum at a level of accuracy sufficient not 
to degrade the overall accuracy of the cosmological pre- 
dictions obtained using the CAMB code, as described in 
Sec. IV Similarly for each tensor fc-mode we set the ini- 
tial conditions at Ne for which y/ii/(l — ei) = 50 and find 
the power spectrum as we evolve forward in time until it 
freezes out at [dlnVh/dlna] < 10^^. Note therefore that 
we do not simply evaluate the power spectrum at hori- 
zon crossing, defined as Cgk = aH, as we discuss below. 
Assuming such an instantaneous freeze-out takes place 
at this time, and further assuming that both tensor and 



scalar modes freeze out nearly simultaneously can have 
notable effects on the estimation of the power spectrum 
variables. 

To first order in slow roll parameters, the equation for 
the scalar power spectrum ( 29 1 becomes [ST) [SS] , 



1 - 2e - 2k + 26 
1 i/2 



X 



87r2M2 cse 



(49) 



Cs k—aH 



where b = 2 — In 2 — 7, and 7 = 0.5772 is the Euler- 
Mascheroni constant. A similar calculation for the ten- 



sor power spectrum (31) is typically evaluated at ten- 
sor mode horizon crossing (fc = aH), however recently 
an approximate expression at scalar horizon crossing 
(csfc = aH) was also given [49l I59j. 



Vh{k) 



[1-2(1-%]^. 



pi 



[l-2(l-6-lnc,)e] 



k=aH 

2H^ 



(50) 
.(51) 



We consider pivot scales, at which spectrum parameters 
are calculated, for scalar and tensor modes as k^,s = 0.01 
Mpc~^ and fc*t = 0.01 Mpc~^, respectively. 

Using the above expressions one can calculate the 
scalar power spectrum normalization, Ag, tilt, n^, and 
running, nrun, of the spectral index, and the tensor spec- 
tral index, rit , at the pivot points, |58[ I59j 



n, - 1 = 



As 

-{2e + r] + k){1 + e + k) 

'f + ^)+2e^ + 2.^,{54) 



(52) 



(53) 



-2b CN 



drig 
d\nk 



(55) 



k=k. 



{2e + ri+ K)(ejv + kjv) + (26^ + 'Hn + kn) 



X {l + e+ k)'^ + 2b(^eNN + 
—2£nn — 2knn, 



2 I 
(56) 



nt 



dlnVh 



(57) 
(58) 



[-2e(l + e + k) + 2(1 - b)eN]k=k,, 
[-2e(l + e + k) + 2(1 - 5 - In Cs)eN + 2eK]k=k,,, 

(59) 



where the approximate expressions in terms of the slow- 
roll parameters are given to second order for Us and nt. 
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and third order for rirun- All parameters are calculated 
at sound horizon crossing, unless stated otherwise. 

The Taylor expanded expressions for the scalar and 
tensor power spectrum, spectral indices and running re- 
duce to previous results for canonical inflation [371 ISSl 
[601 EH and for general inflation [46l EH [Ml [59] to the 
orders quoted in those papers. 

We deflne the tensor-to-scalar ratio without Taylor ex- 
pansion approximations as r exact. 



^exact 



'Ph{k*t)\f 



reeze—out 



l^Q ) I freeze 



(60) 



A common approximation is to calculate r at sound hori- 
zon crossing (cgfc = aH), assuming that the tensor and 
scalar modes freeze out at roughly the same time. To 
second order, this expression is given by [55], 



= 16cse [1 + 2k - h{-q + k) 



(61) 



The approximation above, however, is only valid when 
scalar and tensor modes cross the horizon at similar in- 
stants |51| . Since we keep the speed of sound Cg and its 
dynamical evolution general, we do not assume this apri- 
ori and instead calculate r directly using the ratio T^/i/'P^, 
(60 1, by solving the mode equations for u^ik = fc,s) and 
Vk{k = fc*t), and calculating Vh and Vc^ at freeze-out. 
To flrst order we can write an expression for r^xact as. 



^exact 



[1 - 2eh + 2beh] 
[l-2ec + 25(ec + f 



(62) 



16cse, 



(63) 



where the approximation above means that we have as- 
sumed instantaneous freeze-out of the scalar and ten- 
sor power spectra at their respective horizon crossings. 
Here and Hh are calculated at = aH, and c^, 
ei^, r]Q and are calculated at Csfc*^ = aH. Now for 
Cs(fc*s) < 1, scalar modes leave the horizon at an earlier 
time compared to the tensor modes. So for Cs(fc*s) <C 1 
we expect < Hq, and since {Hh/Ho)'^ is a stronger 
effect than the single factor of e^, we expect therefore 



approx 



< 1. 



Similarly for Cs{k^,s) 3> 1 we expect 
c > 1. We verify these results numer- 
IV A| and flnd that the approximate ex- 



exact I \ 

to get T exact It a 
ically in Sec. 

pression for r can give signiflcant discrepancies from the 
actual tensor to scalar ratio for models in which ^ 1. 

This behavior was recently shown to give a modified 
expression for the tensor-to-scalar ratio ^59J, 



= 16cse [1 + 2k - 6(77 + k) + 2e In c^] 



(64) 



which we find is an excellent analytical approximation 
for our r exact- The fact that r exact differs significantly 
from r approx tells us that, as Cg deviates from 1, the 
next-to-leading order contribution in the expression for r 
becomes important, the Taylor expansion which is often 
used breaks down, and next order terms, as given in (64 1, 
are required. 



IV. OBSERVABLE PREDICTIONS OF 
INFLATION 

We apply the formalism outlined in Sees. [IT] and |III| 
to generate evolutionary trajectories for a general infla- 
tionary model. When constraining flow parameters with 
observational data typically two conditions can be con- 
sidered: 

Condition 1: Constraints on the flow parameters at 
horizon crossing from the form of the observed primordial 
power spectrum. 

Condition 2: The end of inflation arises when e = 1, 
when accelerated expansion as deflned by ([s]) ceases, and 
requires that observable scales crossed the horizon a rea- 
sonable number of e-foldings, say, iVg ^ 50 — 80, before 
the end of inflation. 

Monte Carlo Markov Chain analyses placing con- 
straints on the flow parameters often solely impose Con- 
dition 1, e.g. [371[151[ini[S2, while other analyses addi- 
tionally impose the more restrictive, theoretically moti- 
vated, restriction in Condition 2, e.g. [SHI HI [S3 ESl ED ■ 

In Sec. |IV A[ we consider the properties of inflationary 
evolutionary trajectories, and the resultant power spec- 
tra, under Condition 2. In Sec. |IVB| we apply con- 
straints on the flow parameters from the current WMAP 
5-year data, Sloan Digital Sky Survey (SDSS) Luminous 
Red Galaxies (LRG) galaxy power spectrum and "Union" 
Type la supernovae data sets using Condition 1, and ad- 
ditionally consider the permitted models under the more 
restrictive Condition 2. 



A. Monte Carlo simulations of inflationary 
trajectories 

We randomly select values at the end of inflation (when 
e = 1) of the flow parameters {k, 'A, 'a}, the scale fac- 
tor a(ten) and the sound speed, Cs, within the following 
intervals: 

a{ten) e [5 X 10"^^ 5 X 10^2% 

Cs € [0,2], 

K e [-0.05,0.05], 

i^i^ / e [-0.05,0.05] 



= 



{I < 5), 
{l>5). 



The Hubble constant at this instant is [9J, 
'3.17708 X 10-" 



H{ten) 



a{te 



pi- 



(65) 



(66) 



We also consider models of inflation that allow super- 
luminal propagation of density perturbations, with the 
speed of sound at the end of inflation, Cg e [0, 2]. Faster- 
than-light propagation has been shown to arise in higher 
order QCD corrections, see for example (65] ESI [67], and 
in many other theories, such as [55] [Ml [ZO]- It has 
been shown recently [71 J that superluminal propagation 
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FIG. 1: Monte Carlo sampling of canonical (cs = 1, k = 'a = 0) [upper panels] and general [lower panels] inflation models 
plotted in the (r, ris) (left) and {nrun,ns) (right) planes, truncated at / = 5 in the 'A and 'a flow hierarchies using linear priors 
(black points) and log priors (blue points) on a{ten), Cs, and flow parameters, in the ranges specified in (65 1. For general 
inflation we plot models that have Cs{kfs) G [0, 1]. 



in generic k-essence theories does not lead to the appear- 
ance of closed causal curves (hence they do not violate 
causality) . 

To obtain the scalar perturbation spectrum we evolve 
the parameters back to the time, Ni^i when y^/(l — e — 
K)\N(i — 50, and then evolve Uk forward to freeze-out of 
the power spectrum. The pivot mode crosses the horizon 
at TVe = N^, where k^s = aH/cs\N,- If e [50,80] 
then the observable parameters {As,ns,dns/dlnk} are 
calculated and the trial is recorded. The tensor spectrum, 
and r, are calculated in an analogous manner, by evolving 
back to Nfii, when yh/{^ ~ ^)\Nhi — 50, and evolving 
forward to obtain Vhik^,t). 

In Fig. [T] we constrast the properties of inflationary 
trajectories satisfying Condition 2 in the case of canoni- 
cal inflation (cg = 1, k = 'a = 0), as discussed in [311 172] . 
and in models of general inflation. The introduction of an 
evolving sound speed noticeably alters the distribution of 
spectrum observables arising from the flow trajectories. 



The asymptotic relation ris — 1 « — r/8 that holds for 
canonical inflation is broadened to n^, — 1 « — r/Sc^, and 
the introduction of a non-zero k gives rise to nearly scale- 
invariant models with non-zero running. Allowing super- 
luminal propagation, with Cs > 1, can give scenarios with 
larger tensor-to-scalar ratios [75]. 

For the main analysis in Sec. |IV| we assume linear 
priors on the flow parameters (most consistent with as- 
suming linear priors on the power spectrum observables 
at lowest order). We note, for interest, however that 
Monte Carlo sampling assuming log priors on the flow 
parameters can alter the sampling of allowed models, as 
in the context of canonical inflation [531 El]. The effect 
is more noticeable in general inflationary models where 
introducing log priors can allow larger tensor scenarios 
to be sampled more efficiently. 

In Fig. [2] we demonstrate the difference between the 
exact tensor-to-scalar ratio coming from fully evolving 
both Vi; and Vh, (60 1, and well approximated by (64 1, 
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FIG. 2: The difFerence between r^xact and r approx [left panel] and the dependence on Ca(fc*s) [right panel] for general inflation 
models with Cs{kts) G [0,1] (blue points), and Cs{k^s) G [1,2] (dark red points), for an order 5 Monte-Carlo simulation with 
ranges as given in ( 65 I . 
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FIG. 3: Evolution of Cs [left panel], k [center panel] and Cgk/aH [right panel] for k — 0.01 Mpc~^, for two example trajectories 
with K > (black) and «: < (light blue) at the epoch when current cosmological scales exit the horizon. We find that for all 
viable trajectories, Csk/aH is monotonically decreasing over the course of infiation so that one does not need to be concerned 
about the prospect of multiple horizon crossings, and freezing and thawing of scalar perturbations. 



and the ratio derived by assuming both tensor and scalar 
freeze-out concurrently at scalar horizon crossing (61 1. 
As introduced in Sec. |III C| for sound speeds different 
from 1, and especially as Cg ^0, the discrepancy between 
the two values becomes significant, as much as 50-60%. 
We can understand these deviations by the fact that as 
Cs deviates from 1, tensor and scalar modes leave the 
horizon, and are frozen, at increasingly disparate epochs. 



Since the tilt Ug and running £^re calculated purely 
from the scalar power spectrum (with no reference to 
the tensor power spectrum) , we expect that the values of 
corresponding Us^exact and n™„, exact, found directly from 
the power spectrum, will be similar to the values obtained 
from the numerical expressions (54 1 and (56 1. We have 



verified numerically that the approximate expression for 
rig is as good as the exact calculation to within a few 
percent. The fact that Ug^exact and Ug approx are in good 
agreement also implies that our approximate expression 



for rig to second order is reasonable, and we do not need 
to calcualte Ug to fourth or fifth order. 

As shown in Fig. [3] at the epoch when observable 
modes cross the horizon, k and Cg may be increasing or 
decreasing, Cgk/aH, however, always decreases monoton- 
ically. Therefore once the modes have left the horizon 
(i.e. Cgk/aH < 1) they do not re-enter and we do not 
have to worry about the presence of multiple horizon 
crossings (with the potential for unfreezing and refreezing 
of fluctuations). 



B. Constraints from cosmological observational 
data 



We have included our general inflationary perturba- 
tion code into the CAMB code [75] to evolve background 
equations and first order density perturbations for a fiat 
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universe containing baryons, CDM, radiation, massless 
neutrinos and use CosmoMC [76] to perform a Monte- 
Carlo-Markov-Chain analysis of the model parameter 
space in comparison to current cosmological data. 

In Table |T] we summarize the priors on the flow param- 
eters for five inflationary scenarios we investigate. We 
use linear priors on the flow parameters, {e,K,'A,'a}, up 
to some Imax for 'A and 'a, Cg and ln(10^°As) at horizon 
crossing for fc^^ — 0.01 Mpc^^. is used to calculate 



the value of H{N^,) at horizon crossing using (49 1 and 
(52l, which then gives a{N^,) — Csk^,s/ H\n^. We choose 
truncations e.t I — 2 (CI) and I = 5 (C2) to demon- 
strate the effect of adding in extra degrees of freedom in 
reconstructing the power spectrum. 

Using the approach described in Sec. |TTT] we calcu- 
late the scalar and tensor power spectra for 5 x 10~^ 
Mpc"i < fc < 5 Mpc"^ For the MCMC analysis, we 
purely consider constraints within this range of observ- 
able scales, i.e. Condition 1 of Sec. |IV[ We do not impose 
the stricter requirement of Condition 2, that pertains to 
the full inflationary history, however we do require that 
< e < 1 (and hence that inflation persists) during all 
times from when the initial conditions are set, up to when 
the power spectrum has converged for all observable k- 
modes. We discuss the effect that this additional condi- 
tion has on the parameter constraints below. 

We constrain the models using a combination of cos- 
mological datasets, including measurements of the CMB 
temperature and polarization power spectrum from the 
WMAP 5-year data release |I21 US], the "Union" set of 
supernovae [77], and the matter power spectrum of Lu- 
minous Red Galaxies (LRG) as measured by the Sloan 
Digital Sky Survey (SDSS) [2U[78]. We include the shift 
parameter, asd, to adjust the matter power spectrum as 
discussed in pi]. 

The MCMC convergence diagnostic tests on each sce- 
nario considered are performed on 4 or more chains using 
the Gelman and Rubin "variance of chain mean" / "mean 
of chain variances" R statistic for each parameter. Our 
ID and 2D constraints are obtained after marginalization 
over the remaining "nuisance" parameters, again using 
the programs included in the CosmoMC package. 

In Figs. [4] and [5] we show the ID and 2D marginalized 
posterior probability distributions for the flow parame- 
ters and power spectrum observables, n^, n™„, and r at 
k^s for the canonical and general inflationary scenarios 
studied, and for the "standard" power law inflationary 
model with scale independent running. Comparing the 
constraints on the canonical and general models, we see 
that inclusion of higher order flow parameters can notice- 
ably change the constraints on the power spectrum prop- 
erties. Model C2 significantly opens up the accessible re- 
gion as compared to the standard power law model and 
CI, allowing larger negative running models and large 
tensor amplitudes. That one obtains constraints on all 
higher order flow parameters — separately in C2, 
can be attributed to our truncation of the flow hierarchy 
at / = 5. If one were to further include more higher order 



parameters in the analysis, then this would further open 
up the parameter space for lower-order parameters. 

Note that since the data has only a finite amount of 
information, adding extra degrees of freedom does not 
necessarily lead to a statisticaly important improvement 
in the fit. Table |T] shows a comparison of the fit obtained 
for primordial power spectra based on the flow param- 
eters for different models, with that for the commonly 
assumed canonical power law spectrum with scale inde- 
pendent running of the scalar spectral index. The alter- 
native parameterization of the primordial spectrum does 
not significantly improve the fit with data, the improve- 
ment in does not outweigh the additional degrees of 
freedom added, but rather allows the primordial power 
spectrum to be reconstructed with more freedom. 

Increasing the magnitudes of the higher order flow pa- 
rameters increases the variation of the scalar spectral in- 
dex over the observed scales and significantly boosts or 
diminishes small- and large-scale power to levels incon- 
sistent with CMB and galaxy matter power spectrum 
observations respectively. We show this in Fig. [6] where 
we plot the ID posterior probability distributions for the 
primordial power spectrum for models CI, C2, and G2. 
Specifically, C2 is better able to fit freedom in the power 
spectrum at large and small scales arising from larger 
experimental and cosmic variance errors on those scales. 
We also find that notable degeneracies exist between the 
higher-order flow parameters in C2, reflecting that the 
number of independent degrees of freedom measured by 
the observed primordial power spectrum is less than the 
number of higher-order parameters employed. Future 
small-scale measurements of the power spectrum, for ex- 
ample using Lyman-a observations, will help to reduce 
this interdependency. 

Allowing general inflation, with Cg 1, k ^ 0, as in 
models G1-G3, predominantly alters the tensor ampli- 
tude posterior distribution, consistent with the overlap 
of observational constraints and the Monte Carlo sam- 
pling using linear priors on the flow parameters in gen- 
eral models shown in Fig. [T] The scalar spectral index 
to first order is dependent on 2(— 2e -I- ^A — k) so that 
allowing an evolving sound speed (with k =/= 0) opens up 
the range of ^A that is consistent with observations in 
comparison to canonical models. 

The bounds on the flow parameters that we obtain in 
our analysis arise from two different sets of constraints: 
(i) observations, and (ii) the e < 1 requirement. In canon- 
ical models, we find the constraint on e imposed by the 
e < 1 condition is very similar to that arising from the 
observational constraint, therefore it does not play a ma- 
jor role. This is not true for general models, however. We 
see from ( [61[ ) that the bound on r alone places a rough 
upper bound on the product Cge, but leaves e and in- 
dividually unbounded, as, for example, in [49 . Insisting 
however, that, for consistency, inflation should occur over 
the observable scales in general inflationary models, i.e. 
e < 1 from when the initial conditions for each observable 
mode are set up until when the power spectrum for all 
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Scenario 


Inflation type 




e 


'A 






A(-21nL) 


A(d.o./.) 


CI 


Canonical 


1 


[0,0.5] 


[-0.1,0.1], Ima. 


= 2 








1.14 





C2 


Canonical 


1 


[0,0.5] 


[-0.5,0.5], Ima. 


= 5 








1.18 


3 


Gl 


General 


[0,1] 


[0,0.5] 


[-0.5,0.5], Ima. 


= 2 


[-0.5,0.5] 





1.16 


2 


G2 


General 


[0,2] 


[0,0.5] 


[-0.5,0.5], Iraa^ 


= 2 


[-0.5,0.5] 





0.96 


2 


G3 


General 


[0,2] 


[0,0.5] 


[-0.5,0.5], lrr^a. 


= 2 


[-0.5,0.5] 


[-1.0,1.0], Ima^ = 1 


1.28 


3 



TABLE I: Summary of the parameter ranges investigated for each of the inflationary scenarios in the MCMC analysis in Sec. 
|IVB[ All ranges are for values as the scalar mode fc*3 crosses the horizon. We also show the change in the effective minimum 
X ~ —2 In L, where L is the likelihood, and number of extra degrees of freedom (d.o.f.) in comparison to the fiducial, canonical 
power law primordial power spectrum with scale independent running. 



observable modes 10""' Mpc""'^ ^ ^ ^ 5 Mpc~^ has con- 
verged, introduces a constraint on over the observable 
range. This restricts the value of e{k^s) and over 
and above the observational constraints arising from the 
power spectrum properties at the pivot point. In Fig. [7] 
we demonstrate this by plotting Cs{Ks) vs. e(fc*s) for a 
sampling of flow parameter combinations for model Gl 
that are consistent with observations at the 95% confi- 
dence level. We find that requiring e < 1 induces the 
constraint on and constrains e(A:*s) < 0.03 at the 

95% confidence level. Note that, as seen in Fig.|4j adding 
in an extra flow parameter, ^a, as in model G2, does not 
significantly alter the observational constraints on k in 
the presence of the e < 1 requirement. Extending the 
truncation to include higher nonzero 'A might alter the 
constraints on e, however such additional degrees of free- 
dom are not statistically warranted by the data, so we 
don't consider such models here. 

In Figs. |8] and |9] we show the ID 68% and 95% con- 
fidence levels for the flow parameters and p, X, and 
XCx, for the field choice £x = cj^, as a function of 
the observable comoving mode k{N), that exits the hori- 
zon N e-folds before the end of infiation (the slice at 
fc* = 0.01 Mpc~^ is directly analogous with the ID con- 
straints shown in Fig. |4]), where [15] . 



piN) 

XCx{N) 
X{N) 

Infc(iV) 



I'M 

/9(7V,)exp / 2€{N')dN', 
Jn* 



<N)piN), 
e{N)cs{N)p{N), 



(67) 

(68) 
(69) 



In fc„ - (iV - TV,) - / €{N')dN', (70) 



with 



p{N,) = 



1 - 2e{N,) + 2b (2e(iV,) + r]{N,) 



k{N, 



(71) 



As has been pointed out in previous analyses using power 
law primordial power spectra |79j . the best measured 
modes are around k ~ 0.01 Mpc~^, with large-scale con- 
straints being limited by cosmic variance. However even 



in the scenarios with higher-order fiow parameters al- 
lowed to vary, the observations, in combination with the 
e < 1 condition, impose interesting constraints on the 
fiow parameters across the full range of observable scales. 



C. Reconstruction of viable inflationary 
trajectories 

In this section we consider the viable trajectories that 
satisfy both conditions 1 and 2, namely the spectral prop- 
erties are consistent with observational constraints in the 
k range measurable by CMB and large scale structure 
experiments, and infiation ends with e = 1 around 60 e- 
foldings after observable modes have exited the horizon. 

Fig. [To] shows the results of ~ 150 reconstructed en- 
ergy density evolutions p(0), and kinetic term evolutions 
XCx{(f>), tracing back 100 e-foldings from the end of in- 
fiation, for canonical and general models CI, C2, and Gl, 
in which the power spectrum properties, As,ns,nrun at 
fc*, and the tensor-to-scalar ratio, are consistent with ob- 
servations at the 2(7 level and e = 1 at the end of infiation. 
We set (p = at scalar horizon crossing for fc*, and (j) at 
other epochs is given by 



0(iV) 



N 



^2cs{N')e{N')dN' . (72) 



N, 



We see that introducing higher order fiow parameters 
and/or the possibility of an evolving sound speed open 
up the range of allowed trajectories. Equally, however, 
observational constraints within the single field infiation- 
ary formalism are already starting to tie down the range 
of allowed infiationary histories for both canonical and 
general infiation. 



V. CONCLUSIONS 

The large amount of data available today from CMB 
and large-scale structure surveys can be used to turn 
around the problem of matching infiationary theory to 
observations, to reconstructing the theory from the data 
itself. In this paper we explore the parameter space of 
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FIG. 4: Comparison of the ID marginalized posterior probability distributions for the flow parameters and observables, for 
the cases of power law inflation with scale-independent running (blue full line) and our models CI (black dotted), C2 (black 
dashed), G2 (red dot-dashed) and G3 (red triple dot-dashed). The constraints for model Gl are the same as those for G2, with 
just Cs cutting off at 1, so we do not show them here. 



general single field inflationary models by using a Monte- 
Carlo-Markov-Chain approach in combination with the 
Hubble flow formalism, and constrain the energy density 
and kinetic energy in such a Lagrangian. 

In order to accommodate models that allow the speed 
of sound, Cs, to vary during inflation, it is important to 
take into account the fact that tensor and scalar modes 
cross the horizon at different epochs. This directly affects 



the tensor-to-scalar ratio, which may deviate by as much 
as 50% from the value obtained by assuming that horizon 
crossing epochs are effectively simultaneous. 

We use the full flow parameter evolution equations to 
solve for the scalar and tensor perturbations spectra, and 
subsequently evolve each spectrum until it freezes out, as 
opposed to evaluating it at horizon crossing, in order to 
get precise predictions for the primordial power spectra 
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FIG. 5: The 68% and 95% confidence regions of the (^A,'^A) [left], (r, ris) [center], and (nrun,?is) [right] parameter spaces. 
Constraints for a standard power law spectrum with constant running are shown in dark blue and light blue, as well as 
canonical inflation models CI (black dotted), C2 (black dashed) and general inflation model G2 (red dot-dashed). 
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FIG. 6: Comparison of the ID marginalized posterior probability distributions, showing 68% (dark blue) and 95% (pale blue) 
confidence limits, for the primordial scalar power spectrum for models CI [left], C2 [center] and G2 [right], in comparison to a 
"standard" power law spectrum with scale independent running (black dashed lines). Increasing the number of flow parameters 
increases the freedom with which the spectrum is reconstructed. In particular we can get significantly greater or smaller power 
at small and large scales, least well measured by the CMB and large scale structure data respectively. The bounds on the 
higher-order parameters in model C2 directly arise from constraining this greater or lesser power at the extreme ends of the 
observed scales. 



over the full range of observable scales. We study five dif- 
ferent classes of models of inflation, summarized in Table 
|lj in light of the latest CMB and large-scale structure 
data and show the observational constraints on the flow 
parameters, observed power spectrum, and typically con- 
sidered observables, n^, Urum and r, at the pivot point 
k = 0.01 Mpc-i. 

Including higher-order slow-roll parameters allows a 
higher dimensional flt to the primordial power spectrum, 
with increased power on large and small scales possible in 
comparison to the commonly considered canonical power 
spectrum with scale independent running of the spectral 
index. In models of general inflation, where we allow Cs 
to vary, the condition that inflation should continue (i.e. 
e < 1) on observable scales imposes a natural bound on 



the value of e and k at horizon crossing of the pivot mode, 
and constrains the value of r. 



In the absence of a sound theoretical explanation for 
inflation, the method of action reconstruction holds a lot 
of promise to give valuable directions for the search of 
such a theory. We impose constraints on the energy den- 
sity and kinetic energy terms in the Lagrangian in light 
of current observations. The next step could now be to 
either explicitly construct a general class of allowed La- 
grangians, or study different possible Lagrangians in light 
of these constraints. Once we have a general Lagrangian 
that could consistently explain the observed inflationary 
properties, the elucidation of its theoretical motivation 
should hopefully be that much closer. 
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FIG. 7: While models with large e and small Ca are consistent with observational constraints from the bound on r, imposing the 
condition that inflation persists while the observable scales exit the horizon places an additional restriction on e(A;,s)- The left 
panel shows Cs(fc*s) vs. e{k,s) for models which satisfy observational constraints. Models which additionally satisfy e(fc) < 1 
(red crosses) have an upper limit on the value of e(fc,s), while for larger e(fc«s) (blue squares) the e(fc) < 1 condition is not 
met (these models are rejected in the MCMC analysis). The right panel shows the evolution of e{k) for observable modes, for 
example models which give constraints at the pivot point in agreement with the data, and which satisfy (red, lower two lines) 
or break (blue, upper two curves that reach e = 1) the e(k) < 1 constraint. 
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FIG. 10: Reconstructed energy density p [upper panels] and kinetic term XCx [lower panels] for canonical models CI [left 
panel], C2 [center] and general inflationary model Gl [right] which satisfy both conditions 1 and 2 from Sec. IV Each has 
As, Us, rirun at = fc* and the tensor-to-scalar ratio consistent with observations at the 2a level, and e = 1 at the end of 
inflation with A'^* ~ 50 — 80. The range of represents 100 e-foldings of evolution, with <^ = at scalar horizon crossing for fc, . 
(For canonical models with the usual scalar fleld definition, Cx = c7^, XCx = X.) 



